A self-adaptive mesh method for the Camassa-Holm equation

A self-adaptive mesh method is proposed for the numerical simulations of the Camassa-Holm equation based on its integrable semi-discretization. It is an integrable scheme, possessing the N-soliton solution (see J. Phys. A, 41 355205). Moreover, it is called a self-adaptive mesh method, because the non-uniform mesh is driven and adapted automatically by the solution. Once the non-uniform mesh is evolved, the solution is determined by solving a tridiagonal linear system. Due to these two superior features of the method, the numerical results of the propagation and interactions of soliton and cuspons agree with exact ones very well even by a small number of grid points.